Solution of the Lindblad Equation in the Kraus Representation 
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The so-called Lindblad equation, a typical master equation describing the dissipative quantum 
dynamics, is shown to be solvable for finite-level systems in a compact form without resort to writing 
it down as a set of equations among matrix elements. The solution is then naturally given in an 
operator form, known as the Kraus representation. Following a few simple examples, the general 
applicability of the method is clarified. 

PACS numbers: 03.65.Yz, 05.30.-d, 42.50.Lc 



I. INTRODUCTION 

A quantum system in interaction with an environment, 
that is a larger system with a continuous number of de- 
grees of freedom, displays an irreversible, dissipative, and 
decohering dynamics, which is described by a master 
equation % B, S EL S @, 0, H • The ability of solving the 
master equations plays a considerable role in different 
physical contexts both from practical and fundamental 
points of view. It enables us to discuss the thermalization 
processes and to highlight the problem of individualizing 
the quantum-classical boundary [§]. Master equations 
are also necessary for the realistic description of many 
quantum systems exploited in quantum information pro- 
cessing da HH, E3, m, EI m M, E3. 

The so-called Lindblad equation [18| is the most gen- 
eral form of linear and Markovian master equation pre- 
serving the probability and the positivity of the density 
operator. Although the Lindblad equation can be de- 
rived, under appropriate approximations or in some lim- 
its, from a microscopic model consisting of the quantum 
system and the environment after tracing out the latter 
degrees of freedom 0, 0, H, EL H, @, U , one can even 
start with it as a phenomenological equation 19]. 

Generally speaking, solving master equations requires 
the resolution of a set of differential equations among 
the matrix elements of the density operator with respect 
to a specific basis, usually the eigenstates of the system 
Hamiltonian. Nevertheless, it is generally a cumbersome 
task to list up all of the N 2 equations, N being the di- 
mension of the density matrix (the number of states), 
and if one selects the basis carelessly, many of the equa- 
tions look coupled to each other, even though they are 
not essentially so on the appropriate basis. 

Moreover, even in those lucky cases wherein one can 
find the solution to the set of the differential equations, 



the density operator of the system is given as a collec- 
tion of the matrix elements with respect to the specific 
basis, instead of being expressed in an operator form. It 
is known that the evolution of the density operator is 
given in the Kraus representation [201 ]. but it is not easy 
to construct it from the set of the matrix elements. Over 
the last decade, several efforts have been done in this di- 
rection and several procedures for constructing operator 
solutions to master equations under appropriate condi- 
tions have been proposed [U H3] • 

In this article, we present a procedure for solving 
Markovian master equations directly in an operator form 
in the Kraus representation, instead of giving it by spec- 
ifying its matrix elements. If one works with the matrix 
elements, one needs to write down an N 2 x N 2 matrix 
which couples the A? -2 equations for the matrix elements, 
and then, realizes that many of them are not really cou- 
pled to each other (for relatively simple systems). In the 
approach discussed here, on the other hand, it is not re- 
quired to write such a sparse matrix, where most of its 
elements are zero, but one can realize the minimal set of 
the necessary equations by looking at the relevant part 
of the dissipator of the master equation, associated with 
the "quantum jumps" [Ull. 

The approach is initially sketched in Sec. [II] with two 
simple examples, few-level systems at zero and finite tem- 
peratures, showing the idea and its powerfulness. Then, 
we discuss its general applicability in Sec. [Til] and con- 
clude the article in Sec. IIVI In Appendices [A] and [B] 
we summarize a derivation of master equations, giving 
possible connections between the master equations dis- 
cussed in the main text and microscopic Hamiltonians, 
and present some details about the solution to the second 
example considered in Sec. Ill Bl 



II. SIMPLE EXAMPLES 
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In this section, the resolution of two master equations 
describing the dynamics of two physical systems is pre- 
sented. The first example we discuss refers to a thermal 
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reservoir at zero temperature whereas in the second one 
a finite-temperature reservoir (T ^ 0) is considered. We 
successfully solve the respective master equations build- 
ing up in both cases solutions in the operator form known 
as the Kraus representation instead of expressing them 
in terms of a collection of matrix elements. These two 
examples aim at illustrating the "key ingredients" and 
the basic steps of the general procedure we are going to 
present in Sec. IIIII 

A. A Three-Level System at the Zero Temperature 

Consider a three-level system at the zero temperature 
and suppose its dynamics is described by the following 
master equation 

p{t) = -i[H,p{t)] 

- -r (|{|i><ii,p(*)} - |o)<ijp(t)|i)<oj^ 

- 7 Q{|2)(2|,^)}-|l)(2|p(i)|2)(l|), (2.1) 

where the states \i) (i = 0, 1, 2) diagonalize the Hamilto- 
nian H , 

H= J2 ( 2 - 2 ) 

1=0,1,2 

This master equation can describe the dissipative dynam- 
ics in the triplet sector of two qubits at the zero temper- 
ature, as shown in Appendix IA II In order to solve the 
master equation (|2.ip . observe that it can be rewritten 
in the following form: 

p(t) = Ap(t) + p(t)A* + 7 £ B lP {t)B\ (2.3) 

i=0,l 

with 

A = -i^-l(|l)(l| + |2)(2|), (2.4a) 

B = |0)(1|, B 1 = \l)(2\. (2.4b) 

This rewriting of the master equation naturally suggests 
one to introduce the transformation 

p(t)=e M Pl (t)e Ah . (2.5) 

It is indeed easy to convince oneself that p\{t) evolves 
according to the equation 

Pl (t) = 7 Si(*)ft(*)St(t), Bi(t) = e~ At B t e At . 

i=0,l 

On the other hand, explicitly evaluating the operators 
Bi(t), we obtain 

B (t) = Boe-^- E ^e-^ 2 , B x {t) = B^^-^*, 

(2.7) 



and then (|2.6p is expressed as 

pi(t) = fBopityBle-* + 1 B lPl {t)B{. (2.8) 

It is worth noting that -Bj's are essentially ladder op- 
erators satisfying the following relations, 

B 2 Q = B\ = BxB = 0, BqBi = |0)(2| ^ 0, (2.9) 

so that we get 

B pi(f)Bl = 1 B B lPl {t)B\Bl B lP \{t)B{ = 0. 

(2.10) 

These equations are easily integrated to obtain 

B m (t)B{ = B lPl (0)B\ = B lP (0)B{, (2.11a) 
B oPl {t)B\ = B oP (0)Bl + 7 <S B 1 p(0)B 1 t B ( l. (2.11b) 

It is now straightforward to deduce the expression of 
pi(t), by simply integrating the right-hand side of (|2.8[) 
with l|2.1ip taken into account, 

Pl (t) - p(0) + (1 - e-^)B oP (0)Bl + ytB lP (fl)B\ 

+ (1 - e-T* - -fte-^)B o B lP (0)BlBl (2.12) 

Finally, coming back to the original picture through the 
operator 

e^ = e -iS2i e -7i/2| 2 )(2|+e- iElt e-^/ 2 |l)(l| 

+ e - lBo *|0)(0|, (2.13) 

we get the explicit expression of the density operator at 
an arbitrary time t > 0, 

p(t) = e At p(0)e Att 

+ (1 - e-^)B Q p(fi)Bl + jte-^B lP (0)Bl 

+ (1 - e-* - 1 te-^ t )B Q B 1 p{Q)B\Bl. (2.14) 

We wish to stress that the master equation has been 
solved in the operator form and its solution naturally falls 
into the Kraus representation. As expected, the proce- 
dure presented in this section does not take into account 
at all the set of the equations among the matrix elements, 
which, in this case, involves a 9x 9 matrix (or at least, two 
equations for the diagonal elements and three equations 
for the off-diagonal elements). We only need two coupled 
equations (|2.10|) . Therefore, a great simplification has 
been achieved through the present approach. A key ele- 
ment is to introduce the transformation (|2.5p and rewrite 
the master equation into the form (|2.6[) . The number of 
the necessary equations (two in this case) relies on the 
structure of the right-hand side of (|2.6p . related to the 
so-called "quantum jumps" [23l.[24|. 

In the case under scrutiny, the thermal reservoir is 
at the zero temperature, which results in the vanishing 
derivative in the second equation in (|2.10p . This feature 
makes easier the resolution of the set of the equations in 
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(12.10[) : in fact, one can proceed one by one, solve the sec- 
ond equation first and then the first one after plugging 
the solution of the second into the first. We stress that 
this feature is special for zero-temperature cases and re- 
flects the fact that no transitions into the higher energy 
levels from the lower energy ones occur at the zero tem- 
perature. Similar structures are found in Refs. [22l. |24||. 



B. A Two-Level System at a Finite Temperature 

The second example is a single two-level system at a 
finite temperature, described by the master equation 



p(t) = --Q[a z ,p(t)} 



7+ ( ^{o+°-,p{t)} -<T-p{t)cr+ 



- 7- ( 2 {c7 - <7 +' P(<)} _ a +P^ a - ) ■ ( 2 - 15 ) 

<7z is a spin operator and a± are ladder operators. The 
first term on the right-hand side describes the unitary 
evolution, the second the spontaneous and the induced 
decays from the upper to the lower level, and the third 
the induced transition from the lower to the upper. This 
master equation is derived from a spin-boson model: see 
Appendix IA 21 

As in the previous example, we arrange the master 
equation into 

pi(t) = j + <J^ Pl {t)<j + e-^ + 7_ ( j +( o I (t) ( 7_eT t , (2.16) 
where p\(t) is defined by (|2.5[) with 

A= ' - - 2iQ.)a x , (2.17a) 



7 =7+ +7-, 7 = 7+ 



(2.17b) 



Starting from (|2.16[) and taking into account that a\ = 0, 
it is immediate to reach the following equations: 



a- + pi(t)er_ = 7 + a- + (7_p I (<)cr + o'_e 7 *, 
a-pi(t)a + = 7_cr_(T + /3i(t)(T_o' + e 7t . 



(2.18) 



^a z p{0) P3 (1 - e"^*) + 2ie~^ t/2 sin fit" 
(1 - e-^*) (j±a^p(0)a + + ^a+p(0)a. 



(2.19) 

At this point, it is worth observing that the above ex- 
amples (|2.ip and (|2. 15[) keep the general characteristics of 
the Lindblad master equation, and this suggests the pos- 
sibility of applying the procedure to more general cases, 
which is the subject of the next section. 



III. APPLICATION TO MORE GENERAL 
MASTER EQUATIONS 

So far, we have solved two master equations for finite- 
level systems, one at zero and the other at finite tem- 
perature. Our approach essentially consists in writing 
down a set of equations for suitable operators through 
which it is possible to reconstruct the complete solution 
of the master equation. In other words, by means of 
the resolution of the partial problems given in (|2.10[) and 
(|2.18[) , it is possible to give the complete resolution of 
(|2.ip and (|2.15p . respectively. In this section, we are go- 
ing to analyze a more general situation and identify a set 
of operators which allow us to reproduce the basic steps 
of the approach previously exploited. To better under- 
stand this point, consider the following master equation 



p(t) = -i[H,p(t)) 



1 



m , n £ I 



t 

ran 



(3.1a) 



where 



H = E„\n)(n\ = ^ E n X nn , X„ 



\m) (n| 



In contrast to the previous example, i.e. (|2. 10[) for the 
zero-temperature case, none of the derivatives in (|2.18p is 
vanishing in general. The solution to this set of equations 
is given in (|B8[) . Now, by integrating (|2.16[) with the 
solution (|B8|) taken into account and by going back to 
the original picture, we obtain the solution to the master 
equation (|2.15[) : 

p(t) = ip(0)(l + e-^* + 2e- 7 ' 5 * /2 cosm) 

+ ^a z p(0)a z (1 + e-^' - 2e-^ t / 2 cos fit) 
_ i mat ( _ e -7"*) _ 2ie -^*/2 sinfi ^ Pit) = A P {t) + p(t)rf + ^ lmn X mn p{t)Xl n , (3.2) 

4 \7^ / m.nel 



(3.1b) 

The summations over m and n in the dissipator are 
taken only over the relevant eigenstates to the interac- 
tion with the reservoir, the set of which is symbolically 
denoted by I. This master equation is not the most gen- 
eral master equation of the Lindblad form, but is the 
most general one derived from the microscopic Hamilto- 
nian (|Aip under the Born-Markov approximation in the 
weak-coupl ing, regime [J, 0, H, 0, H, 00, HI or in the van 
Hove limit [25|, [2(| [27j , as far as neither the energy 
nor the energy differences uj mn in I is degenerated (see 
Appendix |A|) . These limitations are just for simplicity, 
and the procedure works for even more general cases (for 
finite-level systems), at least in principle. 

Following the procedure presented in the previous sec- 
tion, rearrange the Liouvillian just in the same form as 
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where 



m,n€:I 

= ~i ^ ( S„ - - Iran J ^"n 
n \ m(El / 



(3.3) 



The transformed operators X mn (t) are easily calculated 
to be 



Y (i\ — e~ At Y e At 

Ytm — 7in)~t ) X„ 

(3.4) 

and the transformed density matrix pi(t) satisfies 

Pl(t) = JmnX mn (t)pi(t)Xl n (t) 
m,n£l 

= lmne^^ lm -' en)t X mn pj(t)Xl n . (3.5) 
From this equation, we get 

£(r m W (e E «i(^-7w)*x mn , Pl (i)X^,) , 



exp ( i(.E m - £ n )i + i y^(7£ 
V eel 



n'el 



where 



(r m )„„' = 8 nn i y^^im — 7£„) + 7„ 



(3.6) 
(3.7) 



and by plugging its solution into (|3.5p . we obtain 

(3.8a) 

which, in the case wherein the inverse T" 1 exists, reduces 
to 



Pi(t) = P(0) 

+ £ lmn[T^{e v ^-l)] nn ,X mn ,p{Q)xl nl . 

(3.8b) 

Since the operator A is already diagonal and therefore 
e M = ^ e -iE n t Xnnj e ^t = J2e^X nn: (3.9a) 



mel 



(3.9b) 



the solution p(i) = e At pi(t)e Ait easily follows. 

From the above procedure for finding the solution to 
the master equation of the Lindblad form, it is clear that 
the remaining problem is to calculate the integral of the 
exponential of the matrix r m defined in (|3 . T[) . Such an 
integral always exists for bounded matrices r m 's. Then, 
the range of effectiveness of this procedure would be ev- 
ident. We understand that this procedure is quite use- 
ful and effective whenever there are only a few channels 
through which the system interacts with the reservoir, 
because in such a case the dimension of the matrix T m , 
which is nothing but the number of levels coupled with 
the reservoir, is small enough to be practically treated, 
irrespectively of the dimension of the system itself, which 
may be very large. On the contrary, as the number of lev- 
els interacting with the reservoir increases, the dimension 
of the matrix T m increases linearly and the procedure 
here expounded becomes less effective. 

The above form of the solution (|3.8[) is very expres- 
sive since it separates the part of the evolution associ- 
ated with the "quantum jumps," X mn > p^X^, , and 



the "no-jump processes," corresponding to e p(0)e y 
[13,[23|. It clearly corresponds to the quantum-jump for- 
mula given in Refs. [ID, [24|], but it should be noted that 
the multiple integrals and the infinite sum of them ap- 
pearing in such a quantum-jump formal formula, which 
are usually difficult to carry out directly, are already 
done in the explicit solution expressed by (|3.8[) . More- 
over, it is worth stressing that the jump contribution is 
evaluated somehow block by block, hence individualiz- 
ing a sort of invariant blocks in the Liouville space. In 
fact, the matrix e Fmt is responsible for the evolution of 
a set of the contractions of the density operator, that is 



-IV 



This is 



(X m opX m0 , X ml pX ml , . . . , X m ^ N _ 1 ^pX m ^ N _ 
well visible for instance from (|3.6p . 

The master equation (|3.ip is derived under the assump- 
tion that there is no degeneracy in the energy difference 
k-W = Em — En" 1 . Possible degeneracies in the energy 
difference would result in a slightly different form of mas- 
ter equation. In such a case, roughly speaking, the oper- 
ator X mn in (|3.1|) has to be replaced with the sum of op- 
erators with the same energy difference and therefore the 
Liouvillian would get additional off-diagonal terms. The 
example (IA14[) bears such terms (the second term of the 
generator). It is worth stressing that the procedure pre- 
sented here equally applies to this case as demonstrated 
in Appendix [B] and further, it works even for the most 
general master equation of the Lindblad form in principle 
(as far as the number of the decay channels is finite). 



A. A Three-Level System at the Zero Temperature 

Finally, it is instructive to consider again the two ex- 
amples solved in Sec. [TTl showing how those solutions are 
reproduced through the general formula (|3.8p . To this 
end, it is enough to find the relevant T m matrices de- 
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fined in (|3.7p . For the first example (|2.1[) , only two decay 
constants are nonvanishing, 701 = 712 = 7, and there are 
only two relevant matrices, T m — (m = 0, 1): 



(3.10) 



They are neither invertible nor diagonalizable. The ex- 
ponentials of these matrices are calculated to be 





1 




- p 


1 









i 


r = 


-1 




















(3.11a) 



(3.11b) 



where (e ro *)n, (e rt)t )i2, and (e Flt )22 elements contribute 
to the solution (|3.8a|) . which reproduces (|2.12p . 



B. A Two-Level System at a Finite Temperature 

For the second example (|2.15|) . we have 
74 



ri 



7 7+ 
7- or 



(3.12) 



v 7- -7 
which are invertible and 

e r ot _ 1_ h-e-'+t + 7+e^ 4 7 +(e 7 -* - e^+ 4 ) \ 

~ l p \7-(eT- f - e-T+*) 7 + e-T+*+7-eT-7 ' 

(3.13a) 

_ _1_ / 7 +e 7+t + 7-e" 7 -* 7+(e 7+t - e^-*) \ 

~ 7^ V7-(e 7 +* - e-T-*) 7_eT+* + 7+e-f-y ' 

(3.13b) 

The elements (e rt)t )i , (e rot )n, (e rit ) 00 , and (e ri ') i are 
relevant in the solution Q3.8p . and we obtain (|2.19p . 



involved in the master equation allows us to find the so- 
lution in a rather simple operator form. In fact, the min- 
imal set of the relevant equations is determined by the 
structure of a part of the dissipator, the generator of the 
evolution of the transformed density operator in (|3.5p . 
which is responsible for the so-called "quantum jumps" 
[23l |24| , and it individualizes a sort of invariant blocks 
in the Liouville space. The solution (|3.8p to the mas- 
ter equation (|3.ip is interesting since it clearly separates 
the "quantum jumps" and the "no-jump processes." A 
remarkable point is that, although we demonstrated the 
present approach with specific examples, one can solve 
general Markovian master equations in principle accord- 
ing to the same idea presented here, as far as the number 
of the decay channels is finite. 
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IV. SUMMARY 

In this article, we have presented a general approach 
to the resolution of the master equations in the Lindblad 
form. Through this procedure, it is possible to solve the 
Markovian master equation in a compact form without 
resort to writing it down as a set of equations among 
the matrix elements. The solution is directly given in 
an operator form, known as the Kraus representation. 
Furthermore, one needs to solve less differential equa- 
tions than for the matrix elements and does not need to 
write down a sparse matrix for the coupled differential 
equations among the matrix elements. This difference is 
remarkable for relatively simple systems like the first ex- 
ample, for which the very simple algebra of the operators 



APPENDIX A: DERIVATION OF THE 
LINDBLAD EQUATIONS 

In this appendix, we wish to give possible relation- 
ships between the master equations discussed in the text 
and microscopic Hamiltonians. Consider a system whose 
Hamiltonian is diagonalized by a set of nondegenerate 
eigenstates, 

H (0) =2 f?(0) |„)(„| =Y J E^X nn , (Ala) 

n n 

and assume that the system interacts with a bosonic 
reservoir, whose annihilation and creation operators are 
denoted by a,k and a k , respectively, and whose Hamilto- 
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man is given by 

H R = J d 3 kuj k ala k , u k > 0. (Alb) 

We consider the general linear interaction between the 
system and the bosonic reservoir, 

flint = ( Xmn 4,n + X L ® Rmn), (Ale) 

^mn = y d 3 kh* mn {k)a k , X mn = \m)(n\, (Aid) 

where the summations over m and n are taken only over 
the eigenstates actually involved in the interaction with 
the reservoir. This set is symbolically denoted by I. 
Note that this interaction Hamiltonian admits counter- 
rotating terms. In this appendix, we further assume that 
there is no degeneracy in the transition frequencies (en- 
ergy differences) 



= - £ (0) 



(A2) 



which requires at least h nn (k) = 0. 

A master equation is then derived from the total 
Hamiltonian 



(A3) 



under the usual Born-Markov approximation in the 
weak-coup ling reg ime P, 0, H, H> IE 2,@L or m the van 
Hove limit [IE Ha S3 » and reads flO) [EEI. The Hamil- 
tonian H is still diagonalized by the same eigenstates as 

H (0) 

even with the Lamb shifts due to the interaction 
with the reservoir taken into account. The decay con- 
stants 7 mn are given by the on-shell form factors with 
the temperature effect of the reservoir included [27j |. 

(A4) 

where 

r m „H = 2n J d 3 k \h mn (k)\ 2 S(uj k - w), (A5a) 

1 



N{uS) 



s /3w _ l ' 



(A5b) 



Note that T mn {oj) = for w < and that the detailed 
balance condition 



In 



(A6) 



holds. 

The master equation (|3.1| is not the most general mas- 
ter equation of the Lindblad form. The characteristic 
feature of (|3.ip is that the superoperators of the unitary 
part (the first term of the generator) H and of the dissipa- 
tor (the second) T> commute with each other, [H, T>] = 0. 
The other types of master equation are derived in differ- 
ent scaling limits [1, [26j . 



1. Two Qubits Immersed in a Bosonic Reservoir at 
the Zero Temperature 



The master equation (|2.ip appears in the dynamics of 
two qubits, say A and B, in interaction with a common 
bosonic reservoir at temperature T = 0. Actually, one 
can consider the total Hamiltonian (|A3|) for the system 
(two qubits) plus the reservoir, with the Hamiltonians 



ffW = ^3+ 5 (4 A V^ +C ri A Vf), (A7a) 
flint = J d 3 k [h*{k)J: + a k + h(k)E-a k ], (A7b) 



and (lAlb|) . where S = er^ + cr^ is (apart fr om the 
normalization) the total spin S. is diagonalized in 

terms of the total spin S (S — 0, 1), 

H^ = QS 3 + g[S(S+l)-S 2 3 -l]= ^fl^l' 



with 



(A8) 

|2> = |T)a|T)b, (A9a) 
|1) = -^(|T)aU)b + U>a|T)b), (A9b) 
|0) = |I)a|I)b, (A9c) 
|s> = ^(|T)aU)b-U)a|T)b), (A9d) 



and 



£<°>=fi, E®=g, E^=-n, E® = -g. 

(A10) 



j 2 — „«, 

The interaction Hamiltonian H\ nt is then written as 



flint = /"dfcv^/i(Jfe)(|0)(l| + |l)<2|)4 + h.c. (All) 

Notice that the singlet sector is decoupled from the reser- 
voir. The triplet sector only suffers from dissipation and 
is described, under the Born-Markov approximation, by 
the master equation (|2.ip . provided the decay constants 
for the two processes, |2) — + |1) and |1) — * |0), given by 

712 = 2T(n-g), 701 = 2T(n + g), (Al2a) 
T(w) = 2tt [ d 3 k \h{k)\ 2 8{uj k - w), (A12b) 



are the same, 712 = 701 = 7- It is understood that the 
energy shifts, say 5i, due to the interaction have already 
been included in Ei of (|2.2p . Actually, we know that 

E 2 = E^ + S 2 , Eh = E{V +6i,E = 4 0) , E s = E^. 

(A13) 



2. Spin-Boson Model at a Finite Temperature 

A master equation of the form 
p(t) = ~lM?z,p{t)\ _ 7o[p(*) - axP{t)<r x ] 



with 



7+ 



~{a+a-,p(t)} - <T-p(t)<r+) 



- 7- { ^{<r-<r+,p(t)} - e+p^o- ) , (A14) 



which reduces to (|2.15[) when 70 = 0, is derived from the 
Hamiltonian (IA3I1 with, for instance 12711. 



# (0) = ^a z , (Al5a) 

Him = <?z J d 3 k[h^(k)a k + h a (k)al] 

+ a x J d 3 k[h* 1 (k)a k +h 1 (k)a{], (A15b) 

and ([Albjl . Note that, in the above general framework 
for deriving (|3.1j) , such terms like the first term of this in- 
teraction Hamiltonian Hi nt , that do not induce transition 
between the two states of the spin system, are excluded, 
i.e. h nn (k) — 0. The second term of the generator of the 
master equation (|A14I) originates from this interaction 
and is not found in (13. ip . The decay constants in the 
generator of (|A14|) and (|2. 15|) are given by [27[ 

70 = /r 1 r o (0+), (Al6a) 

7+ = [i + jv(n)]ri(n), 7 _ = jv(n)ri(n), (Ai6b) 

where 

Ti(w) = 2ir d 3 k\hi(k)\ 2 S(uj k ~uj) (i = 0,l). (A17) 



Between the two decay constants, the detailed balance 
condition holds. 



7 = 



(A18) 



APPENDIX B: SOLVING A MASTER EQUATION 
FOR A TWO-LEVEL SYSTEM 



The master equation (|A14|) for a two-level system at a 
finite temperature is solved as follows. This is the master 
equation (|2 . 15[) when 70 = 0, and is derived from a spin- 
boson model: see Appendix IA 21 

First, we arrange the master equation (|A14[) into 

pi(t) = 10(?zPl{t)<Jz 

+ 7+ CT_pi(£)cr + e~ 7t + 7_er +i oi(i)cr_e 7t 
= (r/oP, + l+e^V- + j-e^V^piit), (Bl) 
where p\(t) is defined by (|2.5[) . i.e. 



p(t) = e M p I (t)e AU = e M Pl (t) 



(B2) 



1 P = 7+ + 7 = 7+ - 7- 



(B3a) 
(B3b) 



and we have introduced superoperators V z , P±, and A, 
defined by 



V z p = a z pa z , V-p = a^pa + , V+p = o+po 

Ap = Ap + pA*. 
These operators TVs satisfy 

vl = 1, v\ = 0, v z r± = v± = r±r z , 

V-V+P- = V-, V+V-V+ = V+. 
Now define p\(t) — e 1aVzt pi(t). then it satisfies 



(B4a) 
(B4b) 



(B5a) 
(B5b) 



-~ Pl (t) = d+e'^V- + ^V+Ypiit). (B6) 

This equation yields the following equations [notice the 
relations (|B5|) ] : 



d_ fe-^V-pi{t)\ _ /-7 7-\ (e^V-piitf 
dt { V-V + pi{t) )~ \ 1+ 0j{ V-V + pi(t) , 

d_ fe^V+pi{t)\ _ / 7 7+ \ fe^V+pi{tY 
dt [v+V-Piit)) J \V+V-p l (t) / 

which are easily solved to give 



(B7a) 
, (B7b) 



-7+* 



yp- 



+ 7_(e 



7-t _ „-1+t 



7+* _ e _ 7-t 



+ 7+(e 7+t 



p(o), 

(B8a) 



p(0). 
(B8b) 



Inserting these results into (|B6[) and integrating its right- 
hand side, we obtain 



ft(*)=P(0) + 4[7+(e 



7-t _ e -7+i 



e 

( 7+e 7-t + 7 _ e -7+t _ y»)p_p + ]p(Q) 



— [7-(- 7+t 



7_e 
P -7-t 



+ ( 7 _eT+ i + 7+e" 7 -' - 7 /3 )7' + ^-]p(0) ! 

(B9) 
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which, together with the relations (|B5[) and e 1oVzt = In order to reach the solution pit) , the following relations 
cosh7oi + V z sinh7 i, yields suffice 

pi(t) = (cosh 70* + V z siDhj t)p(0) 

plot 

+ — [ 7+ (eT-* - e~^)V- 

+ (7+e 7 - 4 + 7 _e-^* - 7 /? )P_P + ]p(0) 

plat 

+ ( 7 _e 7 +* + 7+e" 7 -* - j )V+V-]p(O). 

(BIO) 

I 

e At V- = e -(T-+T°)*P_, e At V+ = e -^+ + ^ t V + (Bll) 

to give 

p(t) = (cosh 7o t + 7> z sinh 7d*)e^*p(0) + -4[7+(l - e~^ Pt )V- + (7+ + 7-e~ 7 "' - 7 /5 e- 7 - t )7>-7>+]p(0) 

'VP 

+ ^[7-(l ~ V+ + (7- + - 7 /3 e- 7 +*)7Wp(0). (B12) 

This reduces to (|2.I9I) when 70 = 0. 
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